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Abstract 



Unitary group branchings appropriate to the calculation of local invariants of density 
matrices of composite quantum systems are formulated using the method of S-function 
plethysms. From this, the generating function for the number of invariants at each degree 
in the density matrix can be computed. For the case of two two-level systems, the gener- 
ating function is F(q) = 1 + q + Aq 2 + Qq 3 + 16g 4 + 23g 5 + 52g 6 + 77q 7 + 150g 8 + 224g 9 + 
396g 10 + 583g n + O (q 12 ). Factorisation of such series leads in principle to the identification 
of an integrity basis of algebraically independent invariants. This note replaces Appendix 
B of our paper[| J Phys A33 (2000) 1895-1914 (quant -ph/0001076) which is incorrect. 
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The measurement problem of detecting nonlocal differences between composite quan- 
tum systems (for example, degrees of entanglement is of great importance for appli- 
cations to quantum computation and communication. At root the question boils down 
to the identification of invariants with respect to unitary transformations which can be 
effected by local operations on each subsystem separately. For the case of two subsystems 
of dimensions Ni and N 2 , the NiN 2 xNiN 2 density matrix p can be regarded in partition 
labelling H as an element of the defining representation {1} of U(Ni 2 N 2 2 ) branching to 
the reducible {!}{l}x{l}{l} representation of U(N 1 )xU(N 2 ) via ({l}x{l})x({l}x{l}) of 
{U{N 1 )xU{Ni))x{U{N 2 )xU(N2)) within {1}x{1} of £7(iVi 2 )x[/(iV 2 2 ). Polynomial invari- 
ants of degree n > are thus U(Ni) xU(N 2 ) singlets of the totally symmetric Kronecker 
power {1} g) {n} = {n} of U (Ni 2 N 2 2 ). According to the standard rules for plethysms[@, U 
the branchings 

U(Ni 2 N 2 2 ) D UiN^) x U(N 2 2 ) D C/(JVi) x U(N 2 ) (1) 
for this plethysm are 

{l}®{n} = ({l}x{l})®{n} = ^Wx{aon}E^WxW 

<t h n a h n 

{^}{^°^} x {A}{Aoa}. (2) 

K \~ Nl n, 
A hjv 2 n 

Here k and A must be iVj., N 2 part partitions of n respectively in order that the cor- 
responding representations of U(Ni) and U(N 2 ) be nonvanishing. However, since the 
product of two representations in a unitary group will contain a singlet if and only if they 
are contragredient, the only singlets occurring are those for which both k o o 3 k and 
A o cr 3 A, or reciprocally a G A o A and a G k o k. The number of singlets F n at degree n 
is thus 

k I-jvi n, 
\\~N 2 n 

where n K \ counts the number of o satisfying this condition"'". The computation thus 
reduces to the evaluation of inner products o of Ni and ^2 part partitions of n (Kronecker 
products in the symmetric group S n ) and leads to the generating function 

00 

F{q) = Y,F n q n . (4) 

n=0 

F(q) is difficult to compute in closed form, but for specific cases can be evaluated to any 
desired degree 0. For example at degree 8 for the 2x2 case, we find 

{6,2}o{6,2} = {8} + {71} + 2{62} + {61 2 } + {53} + 2{521} + {51 3 } + {4 2 } + {431} 

'''Also a should have at most min(iVi 2 , N2 2 ) parts. In Appendix B of our paper jjj], n K \ was erroneously 
taken as 1 



+ {42 2 }; 

{5, 3} o {5, 3} = {8} + {71} + 2{62} + {61 2 } + {53} + 2{521} + {51 3 } + {4 2 } + 2{431} 
+ 2{42 2 } + {421 2 } + {3 2 2} + {3 2 1 2 } + {32 2 1}; (5) 

leading to a contribution (including multiplicity) of ri{6,2},{5,3} = ^{5,3}, {6,2} = 18 to F 8 . In 
this way we calculate 

F(q) = l + g+4g 2 + 6g 3 + 16g 4 + 23g 5 + 52g 6 + 77g 7 + 150g 8 + 224g 9 + 396g 10 + 583g n + O(g 12 ). 

(6) 

This generating function should be compared with the Molien series 0] defined via group 
integration, 



The equivalence between the two series can be readily established in the S"-function for- 
malism. Write g as the element of U(NfN$) corresponding to the adjoint action of 
U(Ni) x U(N 2 ) on p. Characters of group representations are generated by taking traces 
(g n ) of powers of g, and hence are polynomials of the class parameters (x% = exp i<f>i, x 2 = 
expi(f>2, . . .). From the integrand of (0) we have directly 
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where a standard form of the Cauchy product identity has been used[[7| (the complete 
Schur functions S{ n y(x) correspond to one part partitions; for one argument Sr n y(z) = 
z n ). The integrand at degree n thus is indeed the reducible U(N\) x U(N 2 ) character 
corresponding to the n'th symmetrised power {n} of the fundamental representation of 
U(Nf N 2 ), and the invariant integration over ?7(A^i)x?7(iV 2 ) serves to project the identity 
representation^. 

Makhlin[[5| has recently examined the 2x2 case and proposed a concrete set of 18 local 
invariants for mixed state density matrices. Definitive confirmation of the completeness of 
such a set is in principle provided by a factorisation of F(q) which establishes an integrity 
basis presentation of the algebra of invariants in terms of a number of free generators 
together with additional relations. Unfortunately (|[) is not computed to sufficiently 
many terms to deduce an unequivocal factorisation, but for example the form 

(1 + x 4 )(l + x 5 )(l + x 6 ) 4 (l + a; 7 ) 2 (l + a; 8 ) 2 (l + x 9 ) 2 
^ ~ (l-a;)(l-x 2 ) 3 (l-a; 3 ) 2 (l-a; 4 ) 3 
= 1 + x + Ax 2 + Qx 3 + lQx 4 + 23x 5 + 52x 6 + 

+ 11 x 1 + 150a; 8 + 224a; 9 + 398a; 10 + 589a; 11 + 982a; 12 + 0(x 13 ) 

may be noted. This has a denominator set signalling 9 free generators, including 1 at 
degree 1 (the trace of p) and 3 at degree 2 (the traces of the squares of p and of the 

* F(q) in (g) agrees with P(z) quoted by |i| up to degree 11 



reduced density matrices). The total count of 9 is to be expected from the dimensionality 
of the coset manifold SU(A)/S(U(2) x 17(2)), namely 8 = 15 - (8 - 1) plus the overall 
singlet trace. To finite degree it is not possible uniquely to identify the denominator 
factors, and the set of free generators given by G(x) differs from that implied by |J. Of 
course, the saturation of terms in G{x) in (H) compared with F(q) in @ beyond degree 
9 also requires that some of the numerator factors (corresponding to invariant quantities 
whose squares or products are relations in the algebra) should be combined together more 
economically. Nonetheless, the total count in G(x) of 21 invariants (9 denominator plus 
12 numerator quantities), and their degrees, is in agreement with Q. 
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